Recoil Polarization for Delta Excitation in Pion Electroproduction 
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We measured angular distributions of recoil-polarization response functions for neutral pion elec- 
troproduction for W — 1.23 GeV at Q 2 = 1.0 (GeV/c) 2 , obtaining 14 separated response functions 
plus 2 Rosenbluth combinations; of these, 12 have been observed for the first time. Dynamical mod- 
els do not describe quantities governed by imaginary parts of interference products well, indicating 
the need for adjusting magnitudes and phases for nonresonant amplitudes. We performed a nearly 
model-independent multipole analysis and obtained values for ReSi+/Mi+ = —(6.84 ± 0.15)% and 
Re_Ei+ /Mi+ = — (2.91±0.19)% that are distinctly different from those from the traditional Legendre 
analysis based upon Mi+ dominance and £ x < 1 truncation. 

PACS numbers: 14.20.Gk,13.60.Le,13.40.Gp,13.88.+e 



Insight into QCD-inspired models of hadron structure 
can be obtained by studying the properties of the nucleon 
and its low-lying excited states using electromagnetic re- 
actions with modest spacelike four-momentum transfer, 
Q 2 . In the very simplest models, quark-quark interac- 
tions with SU(6) spin-flavor symmetry suggest that the 
dominant configuration for the nucleon consists of three 



quarks in an Estate with orbital and total angular mo- 
menta L = and J = 1/2, while the lowest excited state, 
the A resonance at Ma = 1.232 GeV with J = 3/2, is 
reached by flipping the spin of a single quark and leav- 
ing L = 0. Thus, the pion electroproduction reaction 
for invariant mass W w Ma is dominated by the M\ + 
multipole amplitude. However, the M/^ — Mfq mass split- 
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ting and the nonzero neutron electric form factor clearly 
demonstrate that SU(6) symmetry is broken by color hy- 
pcrfine interactions that introduce D-state admixtures 
with L = 2 into these wave functions lj. Although 
quadrupole configurations cannot be observed directly in 
elastic electron scattering by the nucleon, their presence 
in both wave functions contributes to S\+ and E\+ mul- 
tipole amplitudes for electroexcitation of the A. Addi- 
tional contributions to these smaller multipoles may also 
arise from meson and gluon exchange currents between 
quarks |2| or coupling to the pion cloud outside the quark 
core Recently, it has also become possible to cal- 

culate N to A transition form factors using lattice QCD, 
albeit in quenched approximation ||. 

The relative strength of the quadrupole amplitudes 
is normally quoted in terms of the ratios SMR = 
ReSi+/Mi + and EMR = ReE 1+ /M l+ evaluated for 
isospin 3/2 at W = Ma, but isospin analysis would re- 
quire data for the n7r + channel also. Fortunately, model 
calculations show that the isospin 1/2 contribution to 
these ratios is almost negligible. For example, one ob- 
tains (SMR, EMR) = (-6.71%, -1.62%) for isospin 3/2 
compared with (—6.73%, —1.65%) for the pit channel us- 
ing MAID2003 at Q 2 = 1 (GeV/c) 2 . Therefore, we 
quote results for the pir° channel without making model- 
dependent corrections for the isospin 1/2 contamination. 

Most previous measurements of the quadrupole am- 
plitudes for A electroexcitation fit Legendre coefficients 
to angular distributions of the unpolarized cross section 
for pion production and employ a truncation that as- 
sumes: 1) only partial waves with < 1 contribute and 
2) terms not involving M± + can be omitted. However, 
a more detailed analysis using models shows that nei- 
ther assumption is sufficiently accurate Q. Therefore, 
it is important to obtain data that are complete enough 
for nearly model-independent multipole analysis without 
relying upon sp truncation or Mx+ dominance. 

More detailed information about the nonresonant 
background can be obtained from polarization measure- 
ments that are sensitive to the relative phase between 
resonant and nonresonant amplitudes. This phase infor- 
mation is needed to test dynamical models that attempt 
to distinguish between the intrinsic properties of a res- 
onance and the effects of rescattering. A few previous 
measurements of recoil polarization have been made for 
low Q 2 with the proton parallel to the momentum trans- 
fer [a, 0| j but their kinematic coverage is quite limited. 
Several recent measurements of beam analyzing power 
have also been made [lfl fill ll2T |. Those experiments 
demonstrated that recent dynamical models do not de- 
scribe the nonresonant background well. More gener- 
ally, there are 18 independent response functions for the 
p(e, e'p)ir° reaction, of which half are sensitive to real and 
half to imaginary parts of products of multipole ampli- 
tudes |l3| . In this Letter we report angular distributions 
for 14 separated response functions plus 2 Rosenbluth 



combinations for W = 1.23 GeV at Q 2 = 1.0 (GeV/c) 2 
that are sufficiently complete to perform a phenomeno- 
logical multipole analysis; twelve of these response func- 
tions are obtained here for the first time. Data for a 
wider range of W will be presented later in a more de- 
tailed paper. 

The observables for recoil polarization can be resolved 
into response functions according to 
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where tr is the virtual jN unpolarized cm cross sec- 
tion, A is the beam analyzing power, and P and P' 
are helicity-independent and helicity-dependent polariza- 
tions expressed in terms of longitudinal, normal, and 
transverse basis vectors I = pn, nocgxf, and t oc n x I. 
Here q is the momentum transfer in the lab and pm 
is the final nucleon momentum in the ttN cm frame. 
The response functions, R, depend upon W, Q 2 , and 
cos 9, where 9 is the pion angle relative to q in the cm 
frame; subscripts L and T represent longitudinal and 
transverse polarization states of the virtual photon while 
superscripts include the nucleon polarization component 
and/or a prime for beam polarization, as appropriate. 
Rosenbluth separation of the combinations vtRl+t — 
vlRl+vtRt and vtRl+t = vlR^+vtRt requires vari- 
ation of the beam energy, which was not performed in this 
experiment. The cm phase space is given by i>q = k/qo, 
where k and qo are the pion and equivalent real pho- 
ton momenta, while the kinematical factors vt = 1, 
v T t = e, v' TT = y/1 - e 2 , v L = e s , v LT = y/2e s (l + e), 
v' LT — -y/2es(l — e) are elements of the virtual photon 
density matrix based upon the transverse and scalar (lon- 
gitudinal) polarizations, e = (1 + 2|ptan 2 4^)~ 1 and 
e s = e Q 2 / c L 2 - Finally, (j> is the angle between the scat- 
tering and reaction planes. 

The experiment was performed in Hall A of Jefferson 
Lab using standard equipment described in Ref. [l4| . 
A beam of 4531 ± 1 MeV electrons, with current rang- 
ing between about 40 and 110 /xA, was rastered on a 15 
cm LH2 target. The beam polarization, averaging 72% 
for the first two running periods and 65% for the third, 
was measured nearly continuously using a Compton po- 
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larimeter, with systematic uncertainties estimated to be 
about 1% 

Scattered electrons and protons were detected in two 
high-resolution spectrometers, each equipped with a pair 
of vertical drift chambers for tracking and a pair of scintil- 
lation planes for triggering. Protons were selected using 
the correlation between velocity and energy deposition in 
plastic scintillators and pion production was defined by 
cuts on missing mass and the correlation between miss- 
ing energy and missing momentum. The proton polar- 
ization was analyzed by a focal-plane polarimeter (FPP). 
Detailed descriptions of the FPP and its calibration pro- 
cedures can be found in Refs. 0, 0|. The electron 
spectrometer remained fixed at 14.1° with a central mo- 
mentum of 3.66 GeV/c, while the proton spectrometer 
angle and momentum were adjusted to cover the angular 
distribution. Although the motion of the spectrometers 
was limited to the horizontal plane, the boost from cm 
to lab focuses the reaction into a cone with an open- 
ing angle of only 13° and provides enough out-of-plane 
acceptance to access all of the response functions, even 
those that vanish for coplanar kinematics. Cross sections 
were deduced by comparison with a Monte Carlo model 
of the phase space and acceptance for each setting, in- 
cluding radiative corrections. This model reproduces the 
observed distributions very well |18j . 

The nucleon polarization at the target in the cm frame 
was deduced from the azimuthal distribution for scatter- 
ing in the FPP using the method of maximum likelihood. 
The likelihood function takes the form 

c = n ^(i+z+v-R) (2) 

events 

where £ represents the false (instrumental) asymmetry, 
R is a vector containing the response functions, and r\ 
is a vector of eventwise calculable coefficients that de- 
pend upon kinematical variables, differential cross sec- 
tion, beam polarization and helicity, FPP scattering 
angles and analyzing power, and spin transport ma- 
trix elements. The system of equations derived from 
d\nC/dR m = is solved using an iterative procedure. 
The procedure was tested using pseudodata: a model 
was used to compute response functions for each accepted 
event, the predicted polarization was transported to the 
focal plane using the same transport matrix as for the 
data analysis, and the azimuthal angle in the FPP was 
sampled according to its probability distribution. The 
pseudodata were then analyzed in the same manner as 
real data. We found that the model responses are re- 
covered with fluctuations consistent with the statistical 
uncertainties. We also found that small deviations be- 
tween acceptance-averaged and nominal Q 2 can be com- 
pensated using a dipole form factor. 

Systematic uncertainties due to acceptance normaliza- 
tion, FPP analyzing power, beam polarization, elastic 
subtraction, false asymmetry, and spin rotation matrix 



elements were evaluated by comparing results from re- 
plays differing by a perturbation of the relevant param- 
eter. The propagation of systematic uncertainties for 
fitted Legendre coefficients or multipole amplitudes was 
evaluated using fits to those data sets. Data for W = 1.23 
GeV at Q 2 = 1.0 (GeV/c) 2 are shown in Fig. [T]for bin 
widths of AW = ±0.01 GeV and AQ 2 = ±0.2 (GeV/c) 2 . 
We show Rl+t, Rlt, and Rtt extracted from the <j> de- 
pendence of a with error bars from fitting; the large er- 
ror bars or missing bins for cos 9^0 reflect inadequate 
<fi coverage for this separation, but the phenomenological 
analyses use the actual differential cross sections. The 
bins of cos 9 for polarization were chosen to give approx- 
imately uniform statistics. Inner error bars with endcaps 
show statistical uncertainties and outer error bars with- 
out endcaps include systematic uncertainties. The sys- 
tematic uncertainties in response functions and derived 
quantities are typically small compared with statistical 
or fitting uncertainties. 

Figure ^ a l so shows predictions from several recent 
models: MAID2003 [ICil, DMT SAID 111111: 

and SL H^|. Although the first three response functions 
in column 1 and the last in column 3 have been observed 
before, the other 12 response functions have been ob- 
served here for the first time. The first two columns 
are determined by real parts of interference products 
and tend to be dominated by resonant amplitudes, while 
the last two columns are determined by imaginary parts 
that are more sensitive to nonresonant amplitudes. Thus, 
one finds relatively little variation among models for the 
first two columns and much larger variations for the last 
two, but none provides a uniformly good fit, especially 
to imaginary responses. 

Having factored out leading dependencies on sin 9, the 
response functions should be polynomials in cos 9 of rel- 
atively low order, especially if the assumption of M\+ 
dominance is valid near the A resonance. However, good 
fits over a range of W require additional terms in Rlt, 
Rtt, R'lti R-l+Ti an d Rtt- The green dashed curves fit 
coefficients of Legendre expansions to the data for each 
response function independently, including terms beyond 
Mi 4. dominance as needed; the extra terms have negligi- 
ble effect upon the SMR and EMR values obtained using 
the traditional truncation formulas. Our results for the 
Leg endre analysis are compared with those of Joo et al. 
[25j for CLAS data in the top section of Table U These 
Legendre results for EMR overlap, but our result for SMR 
is more precise and significantly smaller. 

The solid red curves in Fig. ^show a multipole analy- 
sis that varies the real and imaginary parts of all s-wave 
and p-wave amplitudes, except ImAfi_, plus real parts 
of 2— multipolcs. Higher partial waves were determined 
using a baseline model, here based upon Born terms for 
pseudovector coupling. We did not vary ImMi_ because 
all models considered predict that it is negligible for our 
W range, yet experimentally it is strongly correlated with 
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FIG. 1: (Color online) Data for response functions at W — 1.23 GeV and Q 2 = 1.0 (GeV/c) 2 are compared with recent 
models and with fits. The labeling distinguishes L, T, LT, and TT contributions to the unpolarized (0) cross section and to 
transverse (t), normal (n), or longitudinal (1) components of recoil polarization with an h to indicate helicity dependence, if 
any. Linear combinations that cannot be resolved without Rosenbluth separation are identified by L+T. Black dash-dotted, 
dotted, short-dashed, and long-dashed curves represent the MAID2003, DMT, SAID and SL models, respectively. The green 
mid-dashed curves show a Legendre fit while the solid red curves show a multipole fit. 



Im5i_. Note that we could not achieve acceptable mul- 
tipole fits without varying the s-wave amplitudes with 
respect to baseline models and we found that the imagi- 
nary part of 5o+ is especially important. Small improve- 
ments for some of the responses can be obtained by vary- 
ing other ci-wave amplitudes also, but the uncertainties 
in the quadrupole ratios increase because higher partial 
waves are not strongly constrained by these data and cor- 
relations between parameters become more severe. Fits 
starting from the MAID2003, DMT, or SL models are 
practically indistinguishable from those shown, but SAID 
is less suitable as a baseline model because some of its 
i w = 2 amplitudes are too large. The insensitivity of 
quadrupole ratios to the choice of baseline model is shown 
in Table [H Therefore, the multipole analysis provides 
nearly model-independent quadrupole ratios; we choose 
as final the results based upon the Born baseline to min- 
imize residual theoretical bias. 

Both Legendre and multipole analyses reproduce the 
data well but the multipole analysis is more fundamen- 
tal, employs fewer parameters (16 vs. 50), and uses the 
data for all response functions simultaneously while the 
more phenomenological Legendre analysis fits each re- 
sponse function independently and ignores the relation- 
ships between Legendre coefficients required by expan- 
sions of those coefficients in terms of products of mul- 
tipole amplitudes. A more detailed paper is forthcom- 
ing that shows that neither assumption of the traditional 



Legendre analysis (sp truncation and M\+ dominance) 
is sufficiently accurate for data with the present levels of 
completeness and precision. The relative error in the tra- 
ditional Legendre analysis is particularly severe for EMR. 

Recent data on quadrupole ratios for Q 2 < 1.6 
(GeV/c) 2 are compared with representative models in 
Fig. El Note that the MAID2003, DMT, and SL models 
included previous EMR and SMR data in their parame- 
ter optimization. The present result for EMR disagrees 
strongly with the SAID prediction and is nearly identical 
to the data for Q 2 = 0, suggesting that EMR is nearly 
constant over this range. Unlike the somewhat smaller 
CLAS results for EMR, our multipole result does not de- 
pend upon sp truncation or M\+ dominance. Similarly, 
our SMR result is close to those for Q 2 < 0.2 (GeV/c) 2 , 
suggesting that SMR is nearly constant over this range 
also. The stronger Q 2 dependence of lattice QCD calcu- 
lations [5j may arise because the quenched approximation 
misses pionic contributions that are expected to be im- 
portant at low Q 2 . 

In summary, we have measured angular distributions of 
14 separated response functions plus 2 Rosenbluth combi- 
nations for the p(e, e'p)n° reaction at Q 2 = 1.0 (GeV/c) 2 
across the A resonance, of which 12 have been obtained 
for the first time. Dynamical models describe responses 
governed by real parts of interference products relatively 
well, but differ both from each other and from the data 
more strongly for imaginary parts that are more sensi- 
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FIG. 2: (Color online) Red circles, present results Irom mul- 
tipole analysis; open squares, MAMI 0, |26| ; open triangle, 
LEGS filled triangles, MIT cross, ELS A |H; open 
circles, CLAS |25^ . Small horizontal displacements are used to 
reduce clutter. Inner error bars with endcaps show statistical 
and systematic errors; where available, outer error bars with- 
out endcaps include model error. Red, green dashed, blue 
dash-dotted, and cyan dotted curves represent MAID2003, 
DMT, SAID and SL, respectively. Magenta bars show lattice 
QCD results @. 



TABLE I: Quadrupole ratios for Q 2 = 1.0 (GeV/c) 2 . The first 
section reports Legendre analyses and the second multipole 
analyses based upon the specified baseline models. 



method/baseline 


SMR, % 


EMR, % 


xl 


Legendre 


-6.11±0.11 


-1.92 ± 0.14 


1.50 


Joo et al. a 


-7.4 ±0.4 


-1.8 ±0.4 




Born 


-6.84 ±0.15 


-2.91 ± 0.19 


1.65 


MAID2003 


-6.90 ±0.15 


-2.79 ±0.19 


1.67 


DMT 


-6.82 ±0.15 


-2.70 ±0.19 


1.67 


SL 


-6.79 ±0.15 


-2.81 ± 0.19 


1.64 


SAID 


-7.38 ±0.15 


-2.53 ±0.20 


1.85 



"Weighted average of CLAS data for Q 2 = 0.9 (GeV/c) 2 from 



tive to nonresonant mechanisms. None of the theoretical 
models considered provides a uniformly good description 
of the polarization data. We performed a nearly model- 
independent multipole analysis and obtained SMR = 
-(6.84 ± 0.15)% and EMR = -(2.91 ± 0.19)%. The 
traditional Legendre analysis also fits the data well but 
gives distinctly smaller quadrupole ratios, demonstrating 



that its assumptions about the relative magnitudes and 
phases of multipoles are not sufficiently accurate. A more 
detailed presentation of the fitted multipole amplitudes 
will be given in a longer paper. 
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